The muon g-2 discrepancy: errors or new physics? 
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Abstract. After a brief review of the muon g— 2 status, we discuss hypothetical errors in the Standard Model prediction that 
could explain the present discrepancy with the experimental value. None of them looks likely. In particular, an hypothetical 
increase of the hadroproduction cross section in low-energy e+e^ collisions could bridge the muon g— 2 discrepancy, but is 
shown to be unlikely in view of current experimental error estimates. If, nonetheless, this turns out to be the explanation of 
the discrepancy, then the 95% CL upper bound on the Higgs boson mass is reduced to about 130 GeV which, in conjunction 
with the experimental 1 14.4 GeV 95% CL lower bound, leaves a narrow window for the mass of this fundamental particle. 
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INTRODUCTION 

The anomalous magnetic moment of the muon, a^, is one 
of the most interesting observables in particle physics. 
Indeed, as each sector of the Standard Model (SM) con- 
tributes in a significant way to its theoretical prediction, 
the precise a^ measurement by the E821 experiment at 
Brookhaven [1] allows us to test the entire SM and scru- 
tinize viable "new physics" appendages to this theory [2]. 

The SM prediction of the muon g—2 is conve- 
niently split into QED, electroweak (EW) and hadronic 
(leading- and higher-order) contributions: a^^ = 

af-°+af^+af-°+a™°. The QED prediction, computed 
up to four (and estimated at five) loops, currently stands 
at flP = 116584718. 10(16)xlO"''[3], while the EW ef- 
fects provide a™ = 154(2) x 10""[4]. The latest calcu- 
lations of the hadronic leading-order contribution, via the 
hadronic e+e^ annihilation data, are in good agreement: 
= 6909(44) X 10-ii[5], 6894(46) x 10-"[6, 7], 
6921(56) X 10-"[8], and 6944(49) x 10-'^[9]. The 
higher-order hadronic term is further divided into 
two parts: a™° a™°(vp) + fl™°(lbl). The first 

one, -98 (1) X 10"''[6], is the 0{a^) contribution 
of diagrams containing hadronic vacuum polarization 
insertions [10]. The second term, also of 0{a^), is the 
hadronic light-by-light contribution; as it cannot be 
determined from data, its evaluation relies on specific 
models. Recent determinations of this term vary between 
80(40) xlO-'Hll] and 136(25) x 10-^'[12]. The most 
recent one, 110 (40) x 1 0^ " [ 1 3] , lies between them. If we 
add this result to a^°, for example the value of Ref. [6] 
(which also provides the hadronic contribution to the ef- 
fective fine-structure constant, later required for our dis- 
cussion), and the rest of the SM contributions, we obtain 



fl™ = 116591778(61) x lO^'^ The difference with the 
experimental value af = 1 16592080(63) x 10"" [1] 
is Afl^ = a™ _ ^sM ^ +302(88) x lO^^', i.e., 3.4c7 (all 
errors were added in quadrature). Similar discrepan- 
cies are found employing the a^° values reported in 
Refs. [5, 8, 9]. For recent reviews of aj^ see [7, 14, 15]. 

The term 0^''° can alternatively be computed in- 
corporating hadronic T-decay data, related to those of 
hadroproduction in e+e^ collisions via isospin symme- 
try [16, 17]. Unfortunately there is a large difference be- 
tween the e+e^- and T-based determinations of a^}'°, 
even if isospin violation coiTections are taken into ac- 
count [18]. The T-based value is significantly higher, 
leading to a small (^ Iff) Aa^ difference. As the e+e^ 
data are more directly related to the calculation than 
the T ones, all recent analyses do not include the latter 
Also, we note that recently studied additional isospin- 
breaking coiTections somewhat reduce the difference be- 
tween these two sets of data (lowering the T-based deter- 
mination) [19, 20], and a new analysis of the pion form 
factor claims that the t and e+e" data are consistent af- 
ter isospin violation effects and vector meson mixings 
are considered [21]. 

The 3.4a discrepancy between the theoretical predic- 
tion and the experimental value of the muon g—2 can be 
explained in several ways. It could be due, at least in part, 
to an error in the determination of the hadronic light-by- 
light contribution. However, if this were the only cause 
of the discrepancy, fl™°(lbl) would have to move up by 
many standard deviations to fix it - roughly eight, if we 
use the fl™°(lbl) result of Ref. [13] (which includes all 
known uncertainties), and more than ten if the estimate 
of Ref. [12] is employed instead. Although the errors as- 
signed to a™°(lbl) are only educated guesses, this solu- 



tion seems unlikely, at least as the dominant one. 

Another possibility is to explain the discrepancy Aa^ 
via the QED, EW and hadronic higher-order vacuum po- 
larization contributions; this looks very improbable, as 
one can immediately conclude inspecting their values 
and uncertainties reported above. If we assume that the 
g—2 experiment E821 is correct, we are left with two 
options: possible contributions of physics beyond the 
SM, or an erroneous determination of the leading-order 
hadronic contribution (or both). The first of these 
two explanations has been extensively discussed in the 
literature; following Ref. [22] we will study whether the 
second one is realistic or not, and analyze its implications 
for the EW bounds on the mass of the Higgs boson. 



ERRORS IN THE HADRONIC CROSS 
SECTION? 

The hadronic leading-order contribution a^'^ can be 
computed via the dispersion integral [23J 
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where a{s) is the total cross section for e+e~ annihila- 
tion into any hadronic state, with extraneous QED cor- 
rections subtracted off, and s is the squared momentum 
transfer. The well-known kernel function K{s) (see [24]) 
is positive definite, decreases monotonically for increas- 
ing s and, for large s, behaves as m^/(3s) to a good 
approximation. About 90% of the total contribution to 
is accumulated at center-of-mass energies ^ be- 
low 1.8 GeV and roughly three-fourths of is covered 
by the two-pion final state which is dominated by the 
p(770) resonance [17]. Exclusive low-energy e+e^ cross 
sections were measured at colliders in Frascati, Novosi- 
birsk, Orsay, and Stanford, while at higher energies the 
total cross section was determined inclusively. 

Let's now assume that the discrepancy Aa^ = — 
aSM = +302(88) x lO-^\ is due to - and only to - hy- 
pothetical errors in g{s), and let us increase this cross 
section in order to raise a^°, thus reducing Aa^. This 
simple assumption leads to interesting consequences. An 
upward shift of the hadronic cross section also induces 
an increase of the value of the hadronic contribution to 
the effective fine- structure constant at Mz [25], 
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{P stands for Cauchy's principal value). This integral 
is similar to the one we encountered in Eq. (1) for 
a^°. There, however, the weight function in the inte- 
grand gives a stronger weight to low-energy data. Let 



us define a,- = /^'^2 dsfi{s) o{s) (/ = 1,2), where the up- 
per limit of integration is Su < M^, and the kernels are 
Ms)=Kis)/i47:') and/2(*) = [M^J{M^,-s)]/{4aK^). 
The integrals a, with i = 1,2 provide the contributions 
to and Aoi^^{Mz), respectively, from 4m| up to Su 
(see Eqs. (1,2)). An increase of the cross section g{s) 
of the form Aa{s) = ea{s) in the energy range y/s e 
[^/sQ — 5/2, -1-5/2], where e is a positive constant 
and 2m,i + 5/2 < y/s^ < — 5/2, increases ai by 

Afli(V^o,5,e) = ej^^°^gl^2ta{t^)f{t^)dt. If we as- 
sume that the muon g—2 discrepancy is entirely due to 
this increase in (J{s), so that Aai(v^Q,5,e) = Aa^, the 
corresponding increase in Aal^^{Mz) is 



Aa2Wso,5)=Aan " ' . 



(3) 
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The shifts Aa2{y/sQ,5) were studied in Ref. [22] for 
several bin widths 5 and central values ^/sq. 

The present global fit of the LEP Electroweak 
Working Group (EWWG) leads to the Higgs bo- 
son mass Mh = 84+26 GeV and the 95% confi- 
dence level (CL) upper bound M^^ ~ 154 GeV [26]. 
This result is based on the recent preliminary top 
quark mass M, = 172.4(1.2) GeV [27] and the value 
A"had(^z) =0-02758(35) [28]. The LEP direct-search 
95%CL lower bound is = 1 14.4 GeV [29]. Although 
the global EW fit employs a large set of observables, M™ 
is strongly driven by the comparison of the theoretical 
predictions of the W boson mass and the effective EW 
mixing angle sin^0^^' with their precisely measured 
values. Convenient formulae providing the Mw and 
sin^e^^P' SM predictions in terms of Mh, M,, Ao!^l^{Mz), 
and av(Mz), the strong coupling constant at the scale 
Mz, are given in [30]. Combining these two predic- 
tions via a numerical ;|f^-analysis and using the present 
world-average values M„ = 80.399(25) GeV [31], 
sin20^^P' = O.23153(16) [32], = 172.4(1.2) GeV [27], 
as{Mz) = 0.118(2) [33], and the determination 
A«had(^z) = 0.02758(35) [28], we get M„ = 89+^^ GeV 
and M™ = 156 GeV. We see that indeed the M^ values 
obtained from the Mw and sin^O^^' predictions are quite 
close to the results of the global analysis. 

The Mh dependence of a^*^ is too weak to provide Mh 
bounds from the comparison with the measured value. 
On the other hand, Aal^^^{Mz) is one of the key inputs of 
the EW fits. For example, employing the more recent (and 
shghdy higher) value Aa^'^^(Mz) = 0.02768(22) [6] in- 
stead of 0.02758(35) [28], the M« prediction shifts down 
to Mh = 881^ GeV and M™ = 145 GeV. In [22] we 

considered the new values of Aa^^^^{Mz) obtained shift- 
ing 0.02768(22) [6] by Aa2{y/sQ, 5) (including their un- 



certainties), and computed the corresponding new val- 
ues of M™ via the combined j^-analysis based on the 

My, and sm^d^X i^^P^ts (for both Aa<'^^(Mz) and aj^^^ 
we used the values reported in [6]). Our results show 
that an increase Ea{s) of the hadronic cross section (in 
\/s € [v^Q — 5/2, + adjusted to bridge the 

muon g—2 discrepancy Aa^, decreases M™, further re- 
stricting the already narrow allowed region for M^. We 
concluded that these hypothetical shifts conflict with the 
lower Umit M]P when s/sq > 1 .2 GeV, for values of 5 up 
to several hundreds of MeV. In [22] we pointed out that 
there are more complex scenarios where it is possible to 
bridge the Aa^u discrepancy without significantly affect- 
ing M™ , but they are considerably more unlikely than 
those discussed above. 

If T data are incorporated in the calculation of the 
dispersive integrals in Eqs. (1,2), 0"'-° significantly in- 
creases to 7110(58) X 10-"[17], a™°(vp) slightly de- 
creases to — 101(1) X 10~'n6, 15], and the discrepancy 
drops to Aa^ = +89(95) x 10"^', i.e. la. While using 
T data almost solves the Aa^u discrepancy, it increases 
Aa^'^'^(Mz) to 0.02782(16) [34, 17]. In [34] it was shown 
that this increase leads to a low prediction which is 
suggestive of a near conflict with M\p, leaving a narrow 
window for M^. Indeed, with this value of Ao^^|j(Mz) 
and the same above-discussed other inputs of the X^' 
analysis, we find an M™ value of only 133 GeV. 

Recently computed isospin-breaking violations, im- 
provements of the long-distance radiative corrections to 
the decay n^n^Vx [19], and differentiation of the 
neutral and charged p properties [20] reduce to some ex- 
tent the difference between t and e+e^ data, lowering the 
T-based determination of a^°. Moreover, a recent analy- 
sis of the pion form factor below 1 GeV claims that T data 
are consistent with the e+e~ ones after isospin violation 
effects and vector meson mixings are considered [21]. In 
this case one could use the e+e" data below ^ 1 GeV, 
confirmed by the T ones, and assume that Aa^ is accom- 
modated by hypothetical errors occurring above ~ 1 GeV, 
where disagreement persists between these two data sets. 
Reference [22] shows that this assumption would lead to 
M™ values inconsistent with M]p. 

In the above analysis, the hadronic cross section a{s) 
was shifted up by amounts Aa(j) = e(y{s) adjusted to 
bridge Aa^y. Apart from the impUcations for Mf,, these 
shifts may actually be inadmissibly large when compared 
with the quoted experimental uncertainties. Consider the 
parameter e = A(7(5') /o(s). Clearly, its value depends on 
the choice of the energy range [^/sq — 8/2, ^/vq + ^ /2] 
where a{s) is increased and, for fixed it decreases 
when 5 increases. Its minimum value, ~ 4%, occurs if 
c{s) is multiplied by (1 + e) in the whole integration 
region, from Inin to infinity. Such a shift would lead to 
70 GeV, well below M^. Higher values of e are 



obtained for narrower energy bins, particularly if they 
do not include the p-o resonance region. For example, 
a huge e ~ 52% increase is needed to accommodate 
Aa^ with a shift of a{s) in the region from 2m up to 
500 MeV (reducing M^^ to 139 GeV), while an increase 
in a bin of the same size but centered at the p peak 
requires e ~ 8% (lowering M™ to 127 GeV). As the 
quoted experimental uncertainty of a{s) below 1 GeV 
is of the order of a few per cent (or less, in some specific 
energy regions), the possibility to explain Aa^ with these 
shifts Aa(s) appears to be unlikely. Lower values of e 
are obtained if the shifts occur in energy ranges centered 
around the p-co resonances, but also this possibility looks 
unUkely, since it requires variations of a{s) of at least 
~ 6%. If, however, such shifts A(7(i) indeed turn out 
to be the solution of the Aa^^ discrepancy, then M™ is 
reduced to about 130 GeV [22]. 

It is interesting to note that in the scenario where 
Aufi is due to hypothetical errors in (y{s), rather than 
"new physics", the reduced M)^ 130 GeV induces 
some tension with the approximate 95% CL lower bound 
Mh 120 GeV required to ensure vacuum stability 
under the assumption that the SM is valid up to the 
Planck scale [35] (note, however, that this lower bound 
somewhat decreases when the vacuum is allowed to be 
metastable, provided its Ufetime is longer than the age of 
the universe [36]). Thus, one could argue that this tension 
is, on its own, suggestive of physics beyond the SM. 

We remind the reader that the present values of 
sin^0g^' derived from the lep tonic and hadronic ob- 
servables are respectively (sin^e^g"), =0.231 13(21) and 

(sin^e^^P')^ = 0.23222(27) [32]. In Ref. [22] we pointed 
out that the use of either of these values as an input pa- 
rameter leads to inconsistencies in the SM framework 
that already require the presence of "new physics". For 
this reason, we followed the standard practice of employ- 
ing as input the world-average value for sin^Q^^' deter- 
mined in the SM global analysis. Since M™ also depends 
sensitively on , in [22] we provide simple formulae to 
obtain the new values derived from different Mf inputs. 

CONCLUSIONS 

We examined a number of hypothetical errors in the SM 
prediction of the muon g—2 that could be responsible 
for the present discrepancy Aa^ with the experimental 
value. None of them looks likely. In particular, follow- 
ing Ref. [22] we showed how an increase Aa{s) =ea{s) 
of the hadroproduction cross section in low-energy e+e~ 
collisions could bridge Aa^. However, such increases 
lead to reduced upper bounds (lower than 1 14.4 GeV 
- the LEP lower bound - if they occur in energy regions 
centered above 1.2 GeV). Moreover, their amounts 



are generally very large when compared with the quoted 
experimental uncertainties, even if the latter were sig- 
nificantly underestimated. The possibility to bridge the 
muon g—2 discrepancy with shifts of the hadronic cross 
section therefore appears to be unlikely. If, nonetheless, 
this turns out to be the solution, then the 95% CL upper 
bound M™ drops to about 130 GeV. 

If T-decay data are incorporated in the calculation of 
a™, the muon g—2 discrepancy decreases to ~ 1 a. While 

this almost solves Aa^j, it raises the value of Aa!^l^{Mz) 
leading to M™ = 133 GeV, increasing the tension with 
the LEP lower bound. One could also consider a scenario, 
suggested by recent studies, where the t data confirm the 
e+e^ ones below ^ 1 GeV, while a discrepancy between 
them persists at higher energies. If, in this case, Aa^u 
is fixed by hypothetical errors above 1 GeV, where 
the data sets disagree, one also finds values of M™ 
inconsistent with the LEP lower bound. 

It has been suggested [37] that a P-wave electromag- 
netic bound state of tt+tt" (pionium) could enter the dis- 
persion relations through 1% mixing with the p in a way 
that significantly increases the hadronic contribution to 
a^. If so, such a state would give little change to the 
Higgs mass determination and would seem to refute our 
claims. However, Ref. [37] is in error. The required mix- 
ing is actually 0.1, not the erroneous 0.01 claim in [37], 
and such large mixing is not possible. The actual effect 
of pionium on is neghgible [38]. 

If the Afl^ discrepancy is real, it points to "new 
physics", like low-energy supersymmetry where Aa^u is 
reconciled by the additional contributions of supersym- 
metric partners and one expects Mh ^135 GeV for the 
mass of the lightest scalar [39]. If, instead, the deviation 
is caused by an incorrect leading-order hadronic contri- 
bution, it leads to reduced M™ values. This reduction, 
together with the LEP lower bound, leaves a narrow win- 
dow for the mass of this fundamental particle. Interest- 
ingly, it also raises the tension with the Mh lower bound 
derived in the SM from the vacuum stability requirement. 
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